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NOTATIONS 


Mean  wave-applied  torque 


Dimensional  friction  factor 


Drag  coefficient 
Wave  group  velocity 

Darcy  frictional  coefficient 
Nondimensional  local  vorticity 
Nondimensional  local  angle  of  r 
Function  proportional  to  sin 


Function  proportional  t 


Instantaneous  total  drag  force 
Instantaneous  horizontal  net  drag 
Instantaneous  vertical  net  drag 
Instantaneous  horizontal  drag  fore 

Instantaneous  horizontal  drag  fore 
Instantaneous  vertical  drag  force 
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Vorticity  at  the  water  surface 
Mean  vorticity 

Steady  state  mean  vorticity  computed  from 
early  measurements  of  the  angle  of  rotation 
of  the  mill  based  on  the  mean  vorticity 
equation 

Steady  state  mean  vorticity  computed  from 
measurements  of  the  angle  of  rotation  of 
mill  based  on  the  relationship  between 
vorticity  and  angle  of  rotation 

Mean  vorticity  at  t = ta 

Steady  state  mean  vorticity  computed  from 


Sintpartiall 


JE  INSIDE 


The  circulation  cells  characteristic  of  these  bars 
depend  on  the  wave-applied  torque  inside  the  breaker  zone, 
mathematically  represented  by  the  curl  of  the  wave-induced 
force  vector  which  is  a function  of  the  radiation  stress 
tensor  components.  This  study  presents  a comparison  between 
the  total  (integrated  over  the  surf  zone)  wave-applied 
torque  calculated  from  two  independent  data  sources  and  that 
predicted  by  theory.  The  results  indicate  that  the  values 
of  the  total  wave-applied  torque  computed  from  independent 
data  sources  differ  by  a factor  of  two,  but  predicted  values 
are  five  to  ten  times  smaller  than  those  obtained  from 

Values  of  the  dimensional  friction  coefficient 
calculated  from  measurements  of  the  rotation  of  the  water 
column  are  in  agreement  with  the  predicted  values  for  the 
Reynolds  Number  characteristic  of  the  model  and  bottom 


f the  device  used  to  measure  the 
r column  are  examined  theoretically  a 


experimentally. 


CHAPTER  1 
INTRODUCTION 

This  study  addresses  the  problem  of  nearshore  circu- 
lation as  a function  of  incident  wave  and  topographic 
characteristics.  It  is  well  known  that  wave  breaking 
induces  a change  in  the  mean  sea  level  and  eventually 
currents,  as  waves  approach  a shoreline  obliquely. 

The  distribution  of  the  wave  height  and  associated 
set-up  are  responsible  for  the  strength  and  direction  of  the 
currents,  higher  breaking  waves  induce  a larger  set-up 
driving  the  flow  to  areas  of  lower  wave  height  and  smaller 

Qualitative  descriptions  of  wave-induced  nearshore 
circulation  have  been  developed  in  the  last  few  decades; 
however,  theoretical  descriptions  have  become  available  only 
over  the  late  sixties  and  early  seventies.  The  driving 
force  for  nearshore  circulation  is  related  to  the  distri- 
bution of  the  gradients  of  the  radiation  stress  tensor 
components  inside  the  surf  zone  (Longuet-Higgins  and 
Stewart,  1962,  1964;  Bowen,  1967).  Bowen  (1969b)  intro- 
duced one  of  the  first  theoretical  studies  associated  with 
nearshore  circulation  considering  a longshore  variation  of 
the  wave  height  as  the  driving  mechanism  for  the  development 


(Sonu,  1972) , indicated  that  the  streamlines  were  narrow  in 
the  outflow  and  broad  in  the  inflow  probably  due  to  non- 
linear mechanisms.  Also,  Sonu  observed  that  the  circulation 
was  pulsational  and  that  the  mean  surface  slopes  agreed  with 
the  trajectories  of  drogues  placed  in  the  circulation  cells. 
Some  laboratory  observations  (Komar,  1971)  and  numerical 
studies  (Noda,  1974)  were  introduced  later  to  investigate 

Mei  and  Liu  (1977)  presented  a linear  analytical 
approach  to  the  nearshore  circulation  driven  by  periodic 
topographical  effects  along  the  shore,  based  on  Bowen's 
(1969b)  paper.  They  found  that  the  position  of  the  return 
flows  depend  on  the  surf  zone  width  and  on  the  on-offshore 

Several  experimental  studies  have  been  carried  out  to 
investigate  wave  transformation  over  three  dimensional 
bathymetries,  such  as  Whalin  11971)  and  Halker  (1974). 

Mhalin  conducted  a series  of  tests  for  nonbreaking  waves 
over  a bottom  topography  with  parallel  circular  contours. 

He  found  negligible  the  reflection  due  to  the  bathymetry, 
small  the  energy  dissipated  over  the  entire  area  of  measure- 
ment and  large  the  amount  of  energy  involved  in  diffraction 
along  the  wave  crest.  Walker  conducted  a series  of  tests 
for  breaking  and  nonbreaking  waves  over  a bathymetry 
simulating  an  idealized  recreational  surfing  shoal.  He 
found  that  measured  wave  heights  at  the  breaking  point  were 
larger  than  predicted  by  Airy  theory,  also  the  use  of  this 


wave  celerity  produces  a cumulative 
refraction  for  waves  of  finite 


theory  in  prediction 
error  in  predicting  1 

Ebersole  and  Dalrymple  (1979)  developed  and  applied  a 
numerical  model  to  longshore  periodic  features  of  the  type 
investigated  here.  This  model  included  convective 
accelerations,  horizontal  mixing  and  a quadratic  bottom 
friction  term  and  was  applied  to  different  beach  shapes. 

The  authors  concluded  that  the  nonlinear  and  lateral  mixing 
terms  are  important  in  nearshore  circulation  processes, 
particularly  for  irregular  bathymetries. 

This  study  comprises  the  results  of  an  experimental 
investigation  to  examine  transformation  and  breaking  of 
normally  incident  waves  over  one  particular  bathymetry  and 
the  resulting  wave-applied  torque  inside  the  surf  zone.  The 
bathymetry  simulates  the  so-called  transversal  bars  and 
represents  a periodic  form  parallel  to  the  coast  super- 
imposed on  a regular  offshore  sloping  bottom. 

Primary  emphasis  is  focused  on  the  hydrodynamics  of 
the  bar-wave  induced  circulation  rather  than  the  resulting 
sediment  transport;  however,  a few  conclusions  can  be  drawn, 
particularly  with  regard  to  the  maintenance  of  the  bars. 

This  work  highlights  the  wave-induced  force  or  torque 
applied  to  the  water  column  inside  the  surf  zone.  This 
force  is  responsible  for  the  currents  inside  and,  by 
continuity,  outside  the  breakers.  Some  measurements  of 
this  force  have  been  performed  and  compared  with  theory. 


The  measurements  conducted  include  wave  height  distri- 
bution, mean  currents  and  time  dependent  vorticity  in  such 
way  that  the  driving  force  can  be  estimated. 


CHAPTER  2 

GOVERNING  EQUATIONS 

Consider  a Cartesian  system  of  coordinates  as  shown  in 
Figure  (2-1)  where  the  y-axis  is  parallel  to  the  coast,  the 
x-axis  is  directed  seaward  and  the  z-axis  has  its  origin  at 
the  water  surface  pointing  upwards.  The  still  water  shore- 
line is  located  at  x » X]_.  The  Cauchy's  equations  in  the 
x,  y and  z directions  and  continuity  for  an  incompressible 
fluid  are  as  follows: 


where  (u,  v,  w)  are  the  velocity  components,  p is  the 
pressure,  p the  water  density,  g the  gravity  and 


ELEVATION 


components. 


wo-dimensional 


The  direction  o 
f incidence  a wi 
n Figure  2-1.  S 


introduce  her 
irrotationality  o 


number,  k,  approaching  the 
propagation  is  inclined  at  a 
normal  to  the  shoreline  as 

e the  refraction  equation 


the  presence  of  currents  has 
(2-5). 


(2-7) 

neglected  in  Equation 


Depth  Integrated  and  Time  Averaged 


i of  the  difficulties  involved 
in  the  description  of  the  nearshore  hydrodynamics  the  depth 
integrated  and  time  averaged  equations  are  used.  These 
equations  deal  with  mean  quantities  and  reduce  the  problem 
from  three  dimensions  to  two  in  the  horizontal  plane. 


The  time  average  (over  one  wave  period)  removes  the  wave- 
related  fluctuations  in  time.  The  averaging  procedure  and 
further  information  about  these  equations  can  be  found  in 
Phillips  (1966).  The  x and  y momentum  equations  are 

||+u|^+v|£  . "9  |£  + Rx  + Tx  (2_8) 


and  the  depth-integrated  and  time  averaged  continuity 
equation  is 


Here  (U,  V)  are  the  mean  horizontal  velocity  components,  n 
the  mean  water  level,  h the  location  depth  and  (Rx,  Ry)  are 
the  bottom  shear  stress  components.  The  driving  force 
components  (Tx,  T ) are  expressed  in  terms  of  the  radiation 
stress  tensor  components  S^jli  =x,  y;  j =x,  y) 


(2-11) 


the  radiation 


tensor  components 


[ (2n  - 1/2) 


e energy  and  n the  r 
he  celerity  of  the  w 


2 L 1 + sh  (2kh)  1 
= [2  th  (kh)  ]1/2 


(2-16) 

derived  using  linear 


All  of  the  equations  above  1 

The  selection  of  a linear  bottom  shear  stress  model  is 
based  on  the  fact  that  all  data  have  been  collected  under 
low  Reynolds  number  experiments,  and,  in  spite  of  the  fact 
that  the  quadratic  form  provides  a more  realistic  represen- 
tation to  the  actual  stresses,  the  mathematical  complexities 
are  increased  greatly.  In  later  chapters  it  will  be  shown 
that  typical  mean  velocities  and  kinematic  viscosities  are 
about  10  cm/s  and  1 cm  /s,  respectively,  and  considering 
a length  scale  of  100  cm,  a Reynolds  number  of  about 


1 results.  The  linear  representation  f 


where  Aj  is  a (dimensional)  friction  coefficient. 

introduction  of  the  depth  integrated  and  time  averaged 
vorticity  equation  will  be  convenient  for  later  purposes. 

By  Cross-differentiating  with  respect  to  x and  y and 
subtracting  Equations  (2-8)  and  (2-91  yields 


where  w,  the  mean  vorticity  is  defined  as 


Vertical  Vorticity  Distribution 

To  facilitate  later  considerations  of  the  vertical 
distribution  of  the  horizontal  vorticity  component  it  is 
convenient  to  introduce  the  corresponding  governing 
equation.  In  boundary  layer  problems  the  equations  of 
motion  can  be  simplified  by  the  fact  that  the  largest 
changes  in  the  velocity  occur  in  the  vertical  direction. 


XI 

If  Equations  (2-1)  and  (2-2)  are  cross-differentiated  and 
subtracted  and  scaling  arguments  are  considered,  all  the 
gradients  of  the  stresses  in  the  x and  y directions  are 
dropped  when  compared  with  those  gradients  in  the  z direc- 
tion the  equation  of  vorticity  results 

||+u|^+v|2  = v2_“  (2-20) 


where  the  local  vorticity  is  defined  as 


and  v is  the  kinematic  viscosity, 
represented  by 


Comparison  between  the  Applied  Longshore  Thrust 


It  is  well  known  that  as  waves 
the  surf  zone,  breaking  results  in 
energy  introducing  a change  in  the 


the  magnitudes  of  the  gradients  of 
by  Equations  (2-8),  (2-9)  and  (2-11 


^ propagate  obliquely  i 
a reduction  of  the  wav 
radiation  stress  compo 


applied  is  related  to 
those  components  as  shown 


with  no  changes 


the  y-direction 


{-direction  is 


balanced  by  the  change  in  the  mean  water  level, 
absence  of  inertia  terms  and  friction. 


In  the  y-direction,  the  change  of  SX) 
applied  longshore  thrust  on  the  surf 


manifested  as  an 


which  in  the  absence 
the  bottom  friction. 


sr  external  forces  is  balanced  by 
from  Equations  (2-25)  and  (2-17) 


which  was  introduced  by  Longuet-Higgins  (1970)  to  represent 
longshore  currents  in  the  absence  of  lateral  mixing. 

For  a beach  with  a periodic  longshore  variation  the 
forces  induced  by  breaking  waves  are  interpreted  as  an 
applied  torque  and  related  to  the  mean  vorticity  as  shown  in 
Equation  (2-18).  I£  the  inertia  and  external  forces  are 
negligible,  the  torque  is  balanced  by  friction  as  follows: 


This  equation  was  considered  by  Bowen  (1969a)  to  represent 
the  simplest  (linear)  theoretical  approach  to  the  rip 
currents  problem  when  a sinusoidal  longshore  variation  of 
the  wave  height  is  assumed. 


CHAPTER  3 

A small  scale  physical  model  was  designed  and 
constructed  to  study  the  circulation  pattern  and  flow 
properties  over  a simple  bathymetry. 

The  bottom  topography  was  configured  to  represent  an 
alongshore  periodic  bar  formation.  In  order  to  minimize  the 
sidewall  effects  the  model  is  symmetric  about  a line  through 
the  ridge  center  in  the  direction  of  the  wave  propagation  so 
the  wave  energy  converges  toward  the  center  of  the  test 
basin.  The  walls  placed  along  the  centerline  of  the  troughs 
represent  streamlines  for  normal  wave  incidence  and  thus  do 
not  alter  the  flow  pattern.  The  amplitude  of  the  bar  was 
small  so  as  to  be  consistent  with  linear  approximations  and 
it  was  decided  to  construct  the  model  with  a uniform  slope 


Coastal  and  Oceanographic  Engineering  Laboratory  o 
University  of  Florida.  The  model  was  separated  fr 
rest  of  the  basin  by  two  walls  of  0.6  m height.  T 
tion  avoided  the  contamination  of  the  flow  over  th 
the  remainder  of  the  basin. 

The  model  was  14  m long  and  3.12  m wide.  Froi 


shoreside  it  starts  w 


E uniform  slope  a 


followed  by  3.65  m of  wave-like  bathymetry  referred  to  as 
the  “bar  area”  and  finally  an  8.07  m section  with  straight 
and  parallel  bottom  contours  (Figure  3-1).  The  bottom  slope 
is  not  uniform  seaward  of  the  “bar  area,"  but  followed  the 


natural  floor  of  the  basin  from  the  wavemaker  until  a depth 
of  24  cm  is  reached.  Figures  3-la  and  3-lb  show  the  lateral 
and  top  view  of  the  model.  Table  3-1  provides  the  actual 
depth  distribution  and  Figure  3-2  displays  contours  of  equal 


The  bottom  was  composed  of  metal  sheets  firmly  affixed 
to  a wooden  forms  with  the  wave- like  shape.  Minor  irregu- 
larities, of  the  order  of  one-half  cm,  were  noted  after  the 
first  day  of  testing,  probably  due  to  the  swelling  of  the 
wood  due  to  submergence  in  water.  However,  these  irregu- 
larities were  far  offshore  and  no  adjustments  were  made. 
During  the  rest  of  experiments  no  further  changes  were 
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Continued 


NOTE:  (1) 


Where  entries  are  omitted  the  bathymetry  is 
uniform  in  the  y-direction  with  the  depth  in  the 
last  column. 


(2) 


bathymetry  is  symmetric 


Introduction 

The  data  acquisition  program  comprised  a series  c 
tests  in  which  wave  height,  vorticity  and 
measured.  Wave  heights  were  monitored  wit 
and  the  data  recorded  on  chart  paper. 

The  vorticity  data  were  inferred  by  t 
the  angle  of  rotation  produced  by  an  eight  vane  "mill" 
placed  parallel  and  close  to  the  bottom  and  supported  by  a 

a video  camera  above  the  mill.  The  wave-induced  torque  was 
estimated  from  measurements  of  the  time  dependent  vorticity, 
through  the  measurements  of  the  time  dependent  angle  of 
rotation  of  the  mill.  Therefore,  Equation  (2-18)  represents 
the  basis  for  the  measuring  procedure.  The  simplest 
possible  approach  to  estimate  the  torque  assumed  that  the 
nonlinear  terms  were  negligible  so  the  measurements  were 
restricted  to  early  stages  in  time  before  the  nonlinearities 
became  significant.  This  consideration  was  satisfied  by 


limiting 


obtained  by  observation 


ot  floating  disk  trajectories  and  in  a limited  number  of 
locations  by  dye  releases.  The  current  and  vorticity 
information  was  taken  by  a video  camera  placed  above  the 
model  and  recorded  on  tape.  The  camera  and  the  video 
recorder  used  in  the  experiments  were  Sony  Trinicon. 

The  symmetry  of  the  bathymetry  required  measurements 
only  in  one  half  of  the  model. 


Wave  Height  Measurements 

The  wave  height  information  was  obtained  by  f 
meters,  each  supported  by  a metal  tripod, 
were  developed  in  the  laboratory  and  they  are  of  the 
"capacitance  type.”  Briefly,  the  wave  gage  operates  with  an 
insulated  wire  at  a frequency  tunable  in  the  range  of  one  to 
two  megahertz.  The  time  varying  height  of  the  water  column 
creates  a varying  capacitance  which  detunes  the  high 
frequency  oscillator  in  an  approximate  linear  fashion.  The 
signal  is  converted  to  voltage  which  produces  an  analog 
output  proportional  to  the  water  surface  displacement.  Once 
a wire  type  and  length  are  tuned,  no  further  adjustments  are 
required;  however,  calibration  of  each  instrument  was  done 
before  and  after  each  test  and  no  changes  were  observed. 

The  calibrations  were  performed  using  a calibration  tank 
where  the  signal  was  compared  with  a graduated  scale. 

Figure  4-1  displays  the  instrument  and  its  metal  support. 


Figure  4-1.  "Capacitance  Type" 


Gage  and 


centimeter. 


rave  gage  was  recorded  on  chart 
ave  gage  precluded  measurements  of 
depth  was  less  than  one 


Vorticitv  Measurements 


The  vorticity 
dependent  angle  of 
by  a metal  rod  (Figure  4 

the  vanes  were  attached 
same  height  and  material 
and  was  provided  with  tw 
friction  on  the  support  : 
is  a 3.0  by  2.5  cm  metal 


eight  vane  mill  supported 

>f  rigid  and  light  plastic, 
ir  cylindrical  hub  of  the 
lb  is  1.2  cm  in  diameter 
ill  bearings  which  reduce 
base  for  the  support  rod 
which  a plastic  protractor 


was  glued  to  allow  i 
shows  a top  view  of  the  instrument.  In  order  to  avoid 
instability  of  the  device  while  measurementes  were  taken,  a 
weight  was  attached  to  the  top  of  the  rod.  The  top  of  the 
mill  was  located  1.5  cm  above  the  bottom. 

The  motion  of  the  mill  was  recorded  on  tape  and  taken 
by  a camera  placed  1 meter  above  the  instrument  supported  by 
a tripod  as  illustrated  by  Figure  4-4.  The  angle  of  rota- 
tion and  the  time  were  recorded  simultaneously  on  tape.  The 
time  provided  by  a digital  clock  appears  superimposed  with 
the  view  of  the  instruments.  Both  the  angle  of 
and  the  time  were  read  later  on  a television  scr 


rotation 


4-3.  Top  View  of  the  Mill 


Video  Camera,  Tripod 
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Figure  4-5  shows  a view  of  the  screen  during  the  reading 
procedure. 

In  order  to  measure  the  change  of  the  angle  of  rotation 
with  time  in  a particular  location,  the  following  steps  were 
carried  out:  first  the  instrument  and  the  camera  were 

placed  in  the  location  of  interest  and  then  the  wavemaker 
was  turned  on,  the  actual  interval  of  time  during  which  the 
measurements  were  performed  was  from  a minimum  of  2 minutes 
to  a maximum  of  7 minutes.  After  each  measurement,  the 
wavemaker  was  stopped  and  the  circulation  allowed  to  cease 
before  additional  measurements  were  conducted. 


The  velocity  field  was  determined  by  observing  the 
trajectories  of  floating  disks.  The  disks  were  designed  to 
follow  the  mean  currents  and  to  avoid  the  surfing  effect 
inside  the  breaker  zone.  They  were  built  with  a top  layer 
of  plastic  to  provide  strength  and  with  a painted  surface 
for  identification  purposes  and  a lower  layer  of  cork  of 
the  required  thickness  to  keep  the  disk  almost  completely 
submerged.  The  drag  was  increased  by  eight  plastic  vanes 
beneath  the  cork  layer. 


The  diameter  of 
1.4  cm  high.  Figure 


e experiments. 


Television  Screen 


During  Reading 


29 


trajectories  were 


by  a video  camera 


recorded  o 


located  10  m above  the 


CHAPTER  5 

DATA  ANALYSIS  AND  DISCUSSION 

The  purpose  of  this  chapter  is,  through  the  analysis 
and  presentation  of  the  data,  to  provide  insight  into  the 
physics  involved  in  the  problem  of  mean  vorticity  generation 
inside  the  surf  zone.  In  particular,  the  wave-induced 
applied  torque  is  examined.  A method  of  analysis  of  the 
vorticity-related  data  was  developed  and  applied.  Since  the 
vorticity  information  was  taken  at  early  stages  of  wave 
development,  the  wave  stabilization  time  was  measured. 
However,  after  waves  break  the  wave  height  distribution  and 
consequently  the  wave-induced  torque,  inside  the  surf  zone, 
depend  primarily  on  the  bathymetry.  The  wave  height  distri- 
bution and  mean  currents  were  measured  inside  and  outside 
the  breaker  zone.  All  measurements  were  performed  for  three 
different  wave  periods. 


Wave  Height 

Introduction 

This  section  presents  the  results  of  the  measurement 
program  to  determine  the  wave  stabilization  time  (ts)  and 
the  distribution  of  the  steady  state  wave  height. 
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Figure  5-4.  wave  Stabilization  Time  (ts) : 

Period  = 1.9  sj  (a)  Actual  Record  and 
(b)  Wave  Height  Versus  Time;  Upper: 


Measurement  No.  1 ( 


3r:  Measurement  No.  2 


stabilized. 


Figures 
distribution 
considered, 
in  Appendix  A 


unsteadiness  of 


current  pattern. 


ow  the  wave  height 

the  different  wave  periods 

e presented  in  tabular  form 


A common  assumption  in  problems  of  this  type  is  that 
inside  the  surf  zone  the  wave  height  is  proportional  to  the 
local  depth  (H  = yh,  where  y - 0.78).  Figure  5-8  shows  the 
relationship  between  wave  height  and  local  depth  for  the 
measurements  obtained. 


Introduction 


Vorticitv 


This  section  addresses  mean  vorticity  (depth 
integrated)  and  the  downward  diffusion  of  vorticity  from  the 
surface,  inside  the  breaker  zone.  A series  of  measurements 
of  the  time-varying  angle  of  rotation,  8m  (which  is  related 
to  vorticity,  S)  was  performed  inside  the  breaker  zone  using 
the  mill  to  measure  0m  is  described  in  Chapter  4.  The 


the  bottom.  A simple  analytical  model  is  introduced  for  the 
representation  of  mean  vorticity  and  the  applied  torque  and 
in  some  cases  the  friction  coefficients  were  calculated  by 
the  method  of  least  squares.  Also,  a vorticity  diffusion 
model  is  developed  which  represents  the  downward  diffusion 
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Figure  5-5.  Wave  Height  Measurements;  Period  = 1.0  s. 
Actual  (•),  H = h (o)  for  y B 0.78 


Figure  5-6.  Wave  Height  Measurements!  Period  = 1.5  s, 


Figure  5-7.  Wave  Height  Measurements;  Period  - 1.9  s. 
Actual  (•),  H = h (o)  for  y - 0.78 


Figure  5-8.  Measure  Wave  Height  (H)  Versus  Local  Water  (h) 
Depth  Inside  the  Sur£  Zone 
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Type  A is  associated  with  those  cases  where  advection 
was  not  important  and  the  applied  torque  was  balanced  by 
frictional  forces  at  an  early  time.  Type  B represents  the 
case  where  the  local  vorticity  was  increased  by  advection 
from  areas  of  higher  vorticity.  In  Type  C the  local 
vorticity  was  diminished  by  advection  from  areas  of  lower 
level  of  vorticity. 

Mean  Vorticity 

A method  is  introduced  below  to  estimate  the  wave- 
induced  torque.  The  method  is  based  on  Equation  (2-18) 


(I)  (II) 


0]  + a2 


(III)  (IV) 


where  the  forcing  term  A2  can  be  written: 


* ' ~7^L) 


(5-J 


In  order  to  develop  a method  to  determine  the  Have-induced 
applied  torque  (A2)  inside  the  surf  tone  the  following 
assumptions  have  been  made: 

(1)  — The  forcing  term  A2  is  ® function  of  the  radiation 
stress  tensor  and  its  gradients  as  indicated  in  Equation 
(5-2).  For  computation  purposes  A2  was  considered  indepen- 
dent of  time  although  in  reality  it  varies  somewhat  since 
some  unsteadiness  in  the  wave  height  data  was  observed.  A2 
was  estimated  from  early  portions  of  the  measurements  of  the 
time  dependent  angle  of  rotation  of  the  mill  described 

(2)  — The  friction  coefficient  A3  depends  primarily  on 
the  Reynolds  number.  For  computation  purposes  it  was  con- 
sidered as  a constant  and  calculated  when  the  mean  vorticity 
was  approximately  constant  in  the  later  portions  of  the 
data. 

(3)  — Since  it  is  to  be  expected  that  the  convective 
terms  would  be  small  in  the  early  stages  of  vorticity 
development  and  to  reduce  the  complexity  due  to  these  non- 
linear terms,  they  were  omitted  from  Equation  (5-1)  and  in 
the  corresponding  analysis.  This  advantage  is  offset  to 
some  degree  because  of  the  need  to  deal  with  the  unsteady 
phase  of  the  vorticity  equation. 

('ll — The  wave-induced  vorticity  depends  on  the  local 
topographical  characteristics.  For  a beach  form  of  regular 
offshore  sloping  bottom  with  the  superimposed  mild  periodic 
form  parallel  to  the  coast  and  assuming  that  the  bathymetry 


and  the  velocity  field  induced  by  breaking  waves  m 
represented  by 


undulations, 
velocity  due 


amplitude  of  the  t 
vorticity  results 


! functions  of  x,  * is  a function  of  > 
ian  velocities  induced  without  the  bol 
- V are  the  contribution  to  the 
; effect  of  the  undulations  parallel  t 
a small  parameter  proportional  to  the 
>ttom  forms.  As  a consequence  the  mea 


The  first  order  mean  velocities  depend  on  the  wave  angle  of 
incidence.  For  normal  incidence  the  alongshore  component  o 


velocity  (V  1 ')  is  zero  and  the  onshore  me; 
is  also  zero  because  the  onshore  transport 
the  return  flow  satisfying  the 
requirement.  As  a consequence 
mean  vorticity  («/’  ) is  zero, 
angle  of  incidence  than  normal  the  circulation  (order 
affected  by  the  presence  of  longshore  currents.  This 


velocity  (0 
■ balanced  by 


first  order 
different 


11972) . Dropping 


has  been  investigated  numerically  by  Sonu 
the  superscript  and  recognizing  that  Term  III  in  Equation 
(5-1)  can  be  neglected  because 
(a)  the  product 


(b)  the  mean  velocity  V is  particularly  small  at  early 
stages  in  time  since  at  breaking  the  onshore  component 
dominates.  Finally,  the  remaining  mean  vorticity  results 


The  vorticity  can  be  related  to  the  angle  of  rotation 
follows: 


Where  36/ at  represents  the  angular  velocity  of  the  fluid. 
Integrating  Equation  (5-7)  an  expression  for  the  mean 
vorticity  is  obtained 


“ = + <“0  + X~l  B*P  ["  V (t  " ta>^  (5-9) 

introducing  Equation  (5-8)  and  integrating  again 


(t  - ta>] 


“ A1  + 5aJ  (t  ' ta)  + Aj  (“o  ■ 5^’ 
fl  - Exp  [-  -jp  (t  - ta)]j 
:he  following  initial  conditions  were  imposed 


Based  on  Equation  (5-10)  the  values  of  Ai , A-  and  A3 
were  determined  for  each  set  by  the  nonlinear  least  squares 
method  using  the  Newton-Raphson  technique.  The  mean 
quadratic  error  was  defined  as  follows: 


where  6 represents  t) 
of  points  considered, 
following  conditions  i 


measured 


s minimum  when  t 


which  represents  a system  of  three  equations  and  three 
unknowns.  Due  to  the  nonlinearity  of  these  equations  it 
was  necessary  to  obtain  convergence  through  iteration. 


Before  the  calculations  o 
it  is  appropriate  to  mention  s 
the  mill  considered  relevant  a 
responds  to  the  velocity  field 
understood  completely.  The  va 


*2  and  a3  are  described 
e mechanical  properties  of 
this  point.  This  instrument 
i a manner  which  is  not 
3 of  the  mill  may  be  con- 


sidered as  a quadratic  drag  device  and  it  was  found  that 
when  placed  in  an  irrotational  field,  rotation  did  occur. 
Briefly,  the  instrument  senses  gradients  in  the  wave  height, 
which  is  equivalent  to  a gradient  in  the  wave  orbital 
velocities  producing  a net  drag  force  which  induces  a 
rotation.  This  effect  is  particularly  relevant  in  the 
absence  of  breaking  waves.  When  waves  break  the  wave- 
induced  vorticity  overcomes  the  effect  of 
Since  the  measurements  of  the  angle  of 
conducted  from  the  starting 
evaluate  the  wave-induced  applied  torque,  the 
grows  with  time  until  reaching  the  steady  stal 
is  particularly  important  because  t 
do  not  necessarily  break,  in  which 
wave  height  gradient  initiating  th< 
an  extraneous  force  in  the  data.  In  order  to  remove  this 
artificial  contribution,  the  arrival  time  of  the  broken 
wave  (ta)  to  the  location  of  interest  is  necessary.  Further 
details  about  the  mechanics  of  the  mill 
Appendix  D. 

To  evaluate  the  effect 


incoming  waves 
and  introducing 


e presented  i 


f this  artificial  force  c 
3 corresponding  to  station  No.  13 


n Figure  5-10. 


id  the  total  interval  i 
approximately  7s.  Tt 


analysis  \ 
different 

(a|  Friction  w; 
analyzed  considering  that 
coincided  with  the  arrival 


£ time  considered  for  t 
evaluate  the  results  t 


•e  considered: 
neglected  (A3  = 0)  and  th 
iat  the  initiation  of  the 
s first  broken  w 


(bl  Friction  was  neglected  (A3  = 0) 
analyzed  considering  the  actual  arrival  c 


Figure  5-11  shows  the  results  for  Cases  a and  b where 
friction  has  been  omitted,  it  can  be  seen  that  at  a time 
approximately  equal  to  the  wave  stabilization  time  both 
A2-curves  are  essentially  the  same  indicating  that,  for 
reasonable  values  of  ta,  the  artificial  torque  introduced  by 
the  wave-induced  drag  on  the  vanes  of  the  mill,  before  any 
arrives,  does  not  alter  the  estimation  of  A2. 
lat  the  mean  quadratic  error  (e2)  increases  with 
:h  a and  b cases.  This  effect  can  be  associated 
Equation  (5-7),  in  the  absence  of  fric- 
3 friction)  become 
reasonable  value  for  a2 
ems  to  be  about  -10  degrees/s2  because  the  error  is  small 
d the  length  of  time  analyzed  is  large  enough  to  balance 
a effect  of  the  artificial  force. 


tion,  behaves 
significant. 


TIME  (s) 


Angle  of  Rotation,  em,  of  the  Mill  Versus 
Time  Corresponding  to  Run  No.  13; 

Period  = 1.5  s;  Location:  x - 3.20  m, 

y = 0.46  m;  ta  Indicates  the  Arrival  Time 
of  the  First  Broken  Wave  and  ts  is  the 
Wave  Stabilization  Time 


When  friction  is  considered.  Equation  (5-13)  must  be 
solved  simultaneously.  To  illustrate  this  case  consider 
Run  No.  13  again.  As  mentioned  before,  the  nonlinear 
character  of  the  equations  requires  initial  values  for  their 
numerical  solution.  These  values  should  be  as  accurate  as 
possible  so  the  absolute  minimum  of  the  mean  quadratic  error 
can  be  obtained.  To  compute  the  initial  friction  coeffi- 
cient the  later  portions  of  the  data  were  considered.  The 
method  requires  the  absence  of  curvature  in  the  measurements 
of  the  time  dependent  angle  of  rotation  after  few  seconds 
the  motion  was  initiated  and  on  the  previous  known  value  of 
the  forcing  term  a2.  For  the  data  shown  in  Figure  5-10  and 
approximately  constant  mean  vorticity  is  present  after  the 

If  “ 0 (5-14) 


and  from  Equation  5-7  the  friction  coefficient  a2  is  esti- 
mated provided  that  the  steady  state  mean  vorticity  is 
known.  This  value  can  be  calculated  from  the  slope  of  the 
data  presented  in  Figure  5-10,  thus 


(5-15) 


Considering  ( n * h)  ^ h = 3.5  cm,  the  initial  value  of 
the  angle  of  rotation  A^  = 0,  the  initial  torque 
A2  ■ -10  degrees/s  and  30m/3t  = -13  degrees/s  the 


friction  coefficient  ft3  = 2.70  cm/s.  The  simultaneous 
solutions  of  Equations  (5-13)  are  Ax  - 0.14  degrees, 

A2  = -18.86  degrees/s2,  a3  = 2.62  cm/s  and  the  mean 
quadratic  error  7?  = 4.07  degrees2. 

When  this  procedure  was  possible  (data  presented  a 
curvature  near  zero  after  the  first  few  seconds) , a3  was 
calculated;  otherwise,  the  first  few  seconds  were  analyzed 
(ensuring  that  friction  was  not  dominant)  and  only  the 
torque  A2  was  calculated. 

The  instrument  was  placed  about  1.0  cm  from  the  bottom 
for  all  measurements.  Thus,  when  interpreting  the  data  on  a 
depth-averaged  basis  as  has  been  discussed  before  the  calcu- 
lated torque  represents  a fraction  of  that  applied  at  the 

e effect  of  the  distribution  over  depth  of 


the  vorticity  generated  (applied  torque) 
in  a later  section  of  this  chapter. 

The  time  interval  for  which  each  se 
analyzed  depended  on  the  arrival  time  of 
wave  (ta)  and  the  characteristics  of  the 
sets  that  were  well-behaved  the  r 
sidered  ranged  up  to  10.5  s (Run 
the  sets  of  data  were  analyzed  fc 
(Run  No.  157)  because  of  the  complicated  behavior 
data  and  must  be  interpreted  as  an  approximation, 
interval  between  the  initiation  of  the  motion  of  t 
and  the  arrival  of  the  first  broken  wave  (ta)  was 


be  addressed 

: the  first  broken 
time  interval  con- 
ill  time  interval 


ranging  up  to  1.2  s.  Due  to  the  smallness  of  ta,  its  effect 
on  the  calculations  of  A2  and  A3  is  considered  unimportant. 

The  bar  area  was  gridded  into  30  cm  by  30  cm  and  one 
vorticity  measurement  inside  the  surf  zone  was  performed  in 
the  center  of  each  square  in  most  of  the  cases.  The  result- 
ing number  of  stations  is  presented  in  Table  5-3.  Figure 
5-12a,  b , c indicates  the  actual  locations  of  each  station. 


Table  5-3.  Angle  of  Rotation  em  Measurement  Runs 


The  results  of  the  analyses  are  presented  in  Tables 
5-4.  5-5  and  5-6.  Por  those  runs  in  which  values  of  A3,  A; 
and  A3  are  not  included,  the  data  are  presented  in 
Appendix  B but  a meaningful  analysis  was  precluded  due  to 
the  complexity  of  the  data. 

Diffusion  of  Vorticity  Inside  the  Breaker  Zone 

Since  the  development  presented  in  the  preceding 
section  relies  on  depth-integrated  properties  such  as  the 
radiation  stress  tensor  and  its  gradients  as  derived  by 
Longuet-Higgins  and  Stewart  (1964),  the  results  presented 
above  have  been  obtained  under  the  assumption  that  the 


Angle  of  Rotation,  Qm,  of 
Mill  Measurement  Stations 

(c)  Period  = 1.9  s inside 
Figures  the  Run  Numbers  A 
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torque  is  uniformly  distributed  over  depth.  However, 
these  properties  are  truly  depth-dependent  as  presented  by 
Borecki  (1981).  It  is  not  the  objective  of  this  work  to 
develop  the  details  of  the  problem  but  rather  to  justify  the 
assumption  that  within  the  surf  rone  the  torque  is  mostly 
applied  at  the  water  surface  so  values  of  the  Sxx  and  SXy 
components  of  the  radiation  stress  tensor  should  be  calcu- 
lated as  a function  of  the  depth. 

Under  the  considerations  mentioned  above,  the  vorticity 
is  diffused  to  the  bottom  after  a time  interval  inversely 
proportional  to  the  viscosity  in  the  water  column.  Because 
of  the  close  proximity  of  the  instrument  to  the  bottom  the 
torque  calculated  from  the  mean  vorticity  equation  through 
the  measurements  should  be  smaller  than  the  actual  torque, 

idealization  of  the  problem  assuming  that  the  torque  due  to 
the  wave  breaking  is  located  entirely  at  the  surface  so  the 
process  of  diffusion  is  reduced  to  the  boundary  value 
problem  sketched  in  Figure  5-13  where  v represents  the 
viscosity  and  Ti  the  wave-applied  torque. 

The  solution  of  this  problem  can  be  found  through  the 
Laplace  transform  method  and  is 
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5-14.  Nondimensional  Vorticity  Distribution 
Over  Depth  Inside  the  Surf  Zone 


which  represents  a system  of  two  equations  and  two  unknowns 
solved  by  the  Newton-Raphson  Method.  To  illustrate  the 
procedure  consider  Run  No.  13.  The  interval  of  time 
considered  was  7.6  s (Table  5-5) . Since  the  measurements 

utilized  in  Equations  (5-18)  and  (5-19) . The  results  were 


The  results  of  the  analysis  for  those  runs  inside  the 
breaker  zone  are  presented  in  Table  5-7.  Figure  5-15  show 
the  correlation  between  the  mean  applied  torque  A2  and  the 
applied  torque  Ti  for  the  same  runs.  For  the  comparison  t 
applied  torque  A2  was  expressed  in  radians/s  . 


Estimation  of  the  Vorticitv  at  the  Water  Surface 


A limited  number  of  experiments  were  carried  out  in 
order  to  estimate  the  vorticity  at  the  water  surface  and  at 
the  same  time  the  vorticity  close  to  the  bottom  inside  the 


breaker  zone, 
fixed  location 


Two  mills  were  considered,  the  mill  at  a 
was  as  described  previously.  The  other  was 
at  freely  in  the  water  surface.  Figure  5-16a 
area  of  measurement  and  the  relative  location 


recorded  on  video  tape.  The  area  of  measurement  was  about 
45  cm  by  45  cm  placed  inside  the  breakers  and  the  wave 
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5-15.  Correlation  between  the  Mean  Wave- 
Applied  Torque,  A2,  and  the  Wave- 
Applied  Torque  at  the  Water  Surface, 
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Local  values  of  the  mean  wave-induced  torque,  A_,  are 
presented  in  Tables  5-4,  5-5  and  5-6  for  three  different 
wave  periods.  Here,  a numerical  integration  of  A 2 is 
presented.  The  area  of  integration  is  limited  to  the 
zone  of  wave  breaking.  The  method  of  integration  used  was 


he  fourth  order  Simpson's  Rule 
f integral  (I)  of  a function  f 
a,  b)  is  approximated  by 


e specific  interval 


of  partitions  and  x is  the  mean  value  of  x. 

Equation  (5-23)  was  applied  in  x and  y and  the  results 
are  presented  in  Table  5-11  for  the  wave  periods 
investigated.  The  number  of  values  of  A2 , for  the  case 

resulting  value  of  (WT)m  is  considered  as  approximate. 


sr 


=a 


Also,  measurements  for  a particular  disk  were  terminated  if 


the  disk  was  carried  out  of  the  field  of  view  of  the  video 
camera.  Complementary  measurements  of  the  velocity  field 
were  made  through  dye  releases,  to  encompass  the  shallower 


practical.  Dye  measurements  are  complicated  by  dispersion 
effects  and  thus  are  not  of  the  same  quality  as  those 
obtained  with  the  disks.  The  dye  results  are  also  included 
in  Appendix  C. 


Characteristics  of  the  Mean  Current  Pattern 


The  circulation  pattern  for  periods  equal  to  1.5  s and 
1.9  s are  shown  in  Figures  5-20  and  5-21.  The  velocities 
represent  the  average  of  all  measurements  made  at  the 
particular  location.  Those  measurements  indicated  by  dashed 
arrows  correspond  to  dye  releases  and  represent  an  approxi- 
mate value  of  the  mean  velocity  given  the  dispersion  of  the 
dye  particles  inside  the  surf  zone. 


Total  Frictional  Torque  Applied  to  the  Bottom 


The  wave-induced  circulation  pattern  about  the  surf 
zone  can  be  approximate  by  the  following  transport  stream 
function,  y,. 


sin  jX- 
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5-20.  Mean  Current  Measurements:  Circulation 

Pattern;  Period  = 1.5  s;  Floating  Disk 
Trajectories  ( ) and  Dye  Releases  ( — 
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Mean  Current  Measurements:  Circulation 

Pattern;  Period  a 1.9  s;  Floating  Disk 
Trajectories  ( ) and  Dye  Releases  ( ) 
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Nondimensional  Transport  Stream  Function, 


Figure 


where  represents  the  magnitude  of  the  transport  stream 
function,  and  Lx  and  I»„  are  the  dimensions  of  the  rectangle 
which  encloses  the  significant  current  pattern  measured  in 
the  model.  In  addition,  the  no  flow  condition  is  exactly 


satisfied  for  y = 0 , y = Ly,  x = xQ  * Lx  and  is  approxi- 
mately satisfied  for  x = xQ.  This  location  was  chosen 
instead  of  the  shoreline  (x  = x^)  because  computational 
reasons.  The  water  depth  equal  to  zero  introduces  an 
artificial  singularity  of  the  integrated  wave-applied 
torque  due  to  the  simplicity  of  the  representation  of  the 
transport  stream  function  given  by  Equation  (5-24).  The 


approximation  x1  = xQ  does  not  introduce  any  significant 
error  in  the  calculations.  Figure  5-22  shows  the  ratio 
4>/4»o  for  different  values  of  (x  - x0)/Lx  and  y/L„. 

The  transport  stream  function  is  related  to  the 
currents  by  the  following  relationships 
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(5-29) 


Considering  Equation  (5-29),  Equation  (5-30) 


(5-31) 


(5-32) 


Considering  Equation  (5-25) 
following  expressions  are  obtained. 
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Mean  Velocity  Profile  (x-component)  at 
Theory  Based  on^Equation  (5-24)  for 


5-15. 
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where  eg  is  the  steady  state  mean  vorticity  computed  from 
mean  current  measurements,  consequently,  uf  can  be  obtained 
from  the  values  of  presented  in  Table  5-15  and  the 
values  of  Lx  and  L shown  in  Tables  5-13  and  5-14.  Also, 
the  steady  state  mean  vorticity  can  be  calculated  from 
Equation  (5-18),  and  the  data  presented  in  Appendix  B. 

For  illustration  purposes,  consider  Figure  5-10.  since 
30m/3t  = -13.0  degrees/s  the  mean  vorticity  calculated  by 
Equation  (5-18)  is  u ■ -26.0  (degrees/s)  or  e„  = 0.45  1/s. 


CHAPTER  6 

COMPARISON  WITH  THEORY 


Introduction 


In  this  chapter  a comparison  is  made  between  tne 
experimental  and  available  mathematical  representations  o£ 
the  physical  processes,  forming  the  focus  of  this  study. 

In  particular,  an  analytical  representation  of  the  bottom 
bathymetry  and  the  refraction  pattern  is  introduced  con- 
sidering the  wave  field  to  be  represented  by  the  shallow 


water  wave  case.  Finally,  a calculation 
wave-induced  applied  torque  and  the  wave 
are  presented. 


of  the  predicted 
set-up  distribution 


Analytical  Representation  of  the  Bathymetry 
The  bottom  bathymetry  represented  in  the  model  may 


approximated 


1.25,  X1  = .91 


x3  = 6.10  m,  t = 5.19  m,  1 = 2»/Lb  where  Lb  • 3.04  n is  the 
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Theoretical  Wave-Applied  Torque 
The  wave-applied  torque  A2  is  expressed  as 


considering  Equation  (2-11) 
Equation  (6-13)  becomes 


assuming  t 


(6-13) 


Since  the  depth,  sin  o and  wave  height  have  been  expanded  in 
terms  of  the  parameter  |i,  the  terms  of  the  applied  torque 
will  be  examined  in  the  same  way.  The  order  of  the 
variables  is 


(6-15) 


- aJ- 


ay 


therefore,  the  wave-applied  torque  becomes 


5 


(6-23) 


] 


The  individual  terms  in  Equation  (6-14)  (Equations  (6-16) 
to  (6-22) ) are  presented  in  Figures  6-3  to  6-5  and  the  sum 
of  those  terms  given  by  Equation  (6-23)  is  presented  in 
Figure  6-6.  The  term  II  has  been  omitted  in  these  figures 

The  integration  of  Equation  (6-23)  using  Simpson's  Rule 
over  the  breaker  rones  indicated  in  Figure  (5-12  gives  the 
total  wave-applied  torque  predicted  by  theory,  (WT)  , for 
the  three  periods  investigated.  These  results  and  the 
values  of  (WT)m  and  FT  are  included  in  Table  6-2  for 


V Assuming' Breaking  Wave  Conditions; 

Upper:  Distribution  of  IV  Given  by 
Equation  (6-19)  in  1/s*;  Lower:  Distribu- 
tion (6-20)  in  1/s2 
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Distribution  o 


ve-Applied  Torque 


Equation  (6-23) ; Values  Inside  Indicate 
A-  in  1/s2,  Assuming  Breaking  Wave 
Condition 


>"«'  Ry’  a 

i (2-17)  a 


'V  V - 


(6-27) 


(2-17) 


, (6-27)  and  (6-20)  in 
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And  considering  first  the  forcing  term  from 

Equation  (2-12) . the  following  expression  is  obtained 


combining  the 
integrating 


d Equation  (6-30)  a 


(6-35) 


(6-36) 


nD  is  a constant  of  integration  which  can  be  evaluated  at 
the  breaking  point  where  the  wave  height  is  expressed  as 

Hb  ■ T<  "b1’  + hb>  <6-381 


where  the  subscript  (b)  indicates  breaking  conditions. 
Substituting  Equation  (6-37)  in  Equation  (6-36)  the  constant 


Finally,  Equation  (6-36) 


which  represents  the  well  known  set-up  distribution  over  a 
sloping  beach. 

Considering  first  the  forcing  term  in  the  y-direction 
and  recalling  Equations  (2-11)  and  (2-12),  it  follows 


and  from  Equation  (6-10) 

,s(!)  , 

-5x“  = M •*»  JS  t (nU>  + h(l,)f3]  (6-43) 

Substituting  Equation  (6-41),  considering  Equations  (6-42) 
and  (6-43) , recalling  Equation  (6-10)  and  operating, 

Equation  (6-32)  can  be  written  as 


(6-45) 


16(3“'  + h1 

integrating  Equation  (6-44 


»gU  + V» 


where  B results  to  be  zero  or  a constant  as  will  be 
demonstrated.  As  Equation  (6-46)  is  substituted  in 
Equation  (6-31)  and  integrated  between  y = 0 and  y = Lv. 

In  order  to  illustrate  this  result.  Figure  6-7  shows 
the  first  order  set  up,  n* , and  the  total  set-up, 
n*1'  + V n*2*  , for  an  arbitrary  location  of  the  breaker  line 


following  parameters 


Wave  Set-Up  Distribution:  A Qualita- 
tive Description;  the  Superscript  (i) 
Indicates  First  Order  Set-Up,  Only  a 
Function  of  x;  the  Superscript  (2) 
Indicates  the  Contribution  Due  to  the 
Bottom  Undulations;  Values  of  the 
Set-Up  in  m. 


3.65 


1X5 


Computed  and  Predicted  Friction  Coefficient 

The  average  value  of  the  friction  coefficient  A3 
resulted  equal  to  1.62  cm/s.  In  order  to  evaluate  the 
representativity  of  that  value  consider  Laminar  flow  which 
is  appropriate  for  the  conditions  given  in  the  model.  The 
linearization  of  the  bottom  shear  stresses  presented  in 
Chapter  2 implies  that 


where  f is  the  Darcy  friction  factor  and  the  wave  orbital 
velocity.  For  laminar  flow  f can  be  expressed  as 


where  Re  represents  the  Reynolds  number  given  by 


Re  = (6-49) 

where  u/o  represents  the  excursion  of  the  water  particle. 
Substituting  Equation  (6-49)  in  Equation  (6-48)  the  friction 
coefficient.  A3,  becomes 
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Considering  the  kinematic  visocity  ranging  from 
to  1 cn)2/s  and  assuming  a wave  period  equal  to  1 


which  is  within  the  range  of  the  computed  values 
presented  in  Tables  5-5  and  5-6. 


CHAPTER  7 

CONCLUSIONS  AND  RECOMMENDATIONS 


A number  of  conclusions  and  recommendations  are 
warranted  as  a result  of  this  investigation.  Conclusions 
drawn  from  this  study  must  be  interpreted  with  a view  toward 
the  effects  of  a more  general  type  of  bottom  bathymetry  in 
spite  of  the  fact  that  the  experimental  data  were  acquired 
for  only  one  topography.  The  conclusions  and  recommenda- 
tions of  this  investigation  are  enumerated  below. 

The  general  circulation  pattern  is  in  qualitative 

the  radiation  stress  tensor.  This  includes  a shoreward- 
directed  flow  over  the  bar  crest  and  a return  flow  from  the 
trough  regions.  The  theoretical  results  presented  in  this 
study  and  those  obtained  by  Mei  and  Liu  (1977)  predicted  a 
counter-circulation  within  a small  region  near  the  shoreline 
which  was  not  observed  in  the  model  data.  This  effect  might 
be  due  to  the  absence  of  lateral  mixing  in 


considerations. 

Measured  values  of  the  wave  height  are  in  good  agree- 
ment with  prediction,  based  on  the  simple  spilling  breaking 
assumption.  The  proportionality  factor  y between  the  wave 
height  and  the  local  water  depth  varies  from  .6  to  1.0.  The 
smaller  values  corresponded  to  regions  of  mild  bottom  slope 
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and  the  larger  values  corresponded  to  regions  of  steep 
bottom  slope , which  is  in  agreement  with  the  results 
obtained  by  Dally  et  al.  (1985). 

The  local  wave-induced  torque  computed  from  measure- 
ments of  the  angle  of  rotation  of  the  mill  indicate  an 
overall  mean  value  of  A2  - -.15  1/s  . Maximum  values  of  A2 
were  obtained  between  the  crest  (centerline)  of  the  "bar" 
and  the  midpoint  between  the  crest  and  trough.  Theory 
predicts  maximum  values  of  a2  at  the  midpoint  between  crest 
and  trough. 

Values  of  the  friction  (dimensional)  coefficient,  A3, 
from  measurements  of  the  angle  of  rotation  of  the  mill, 
yield  an  average  value  of  A3  - 1.5  cm/s,  which  is  within  the 
expected  range  based  on  Reynolds  Number  characteristic  of 


Based  on  the  transient  vorticity  diffusion  model  the 
wave-applied  torque  T j/h  located  at  the  water  surface  is 
about  two  times  larger  than  the  wave-applied  torque,  A2, 
considered  to  be  uniformly  distributed  over  depth. 

Measurements  (steady  state)  of  the  angle  of  rotation  of 
the  mill  at  the  water  surface  (“free"  mill)  and  at  mid-depth 
("fixed"  mill)  indicate  that  the  vorticity  at  the  water 
surface  is  about  two  times  larger  than  that  at  mid-depth, 
which  is  consistent  with  that  mentioned  above. 

The  total  (integrated  over  the  surf  zone)  wave-induced 
torque  was  computed  from  two  independent  data  sources  and 
compared  with  the  theory  of  Mei  and  Liu. 


The  total  (transient)  wave-induced  t 

the  angle  of  rotation  o 
integrated  and  time  averaged,  over  one  wave  period! 
vorticity  equation  in  the  absence  of  nonlinear  terms. 

The  total  (steady  state!  frictional  torque  applied  to 
the  bottom,  FT,  was  comouted  from  mean  current  measurements 
based  on  the  mean  vorticity  equation  and  assuming  an 
aporoKimate  form  of  the  transport  stream  function  in  accord 
with  observations  of  the  circulation  pattern. 

Since  (WT!m  and  FT  should  be  equal  from  theory,  the 
resulting  values  were  compared  along  with  the  predicted 
total  wave-induced  torque  (V!T)_,  for  the  three  periods 
investigated.  The  results  indicate  that  ( WT ) m and  FT 
(which  are  calculated  from  independent  data  sources!  are  in 
reasonable  agreement.  However,  the  predicted  values  of 
!WT)p  are  smaller  by  a factor  of  five  to  ten  than  (WT!m  and 
FT.  Ho  explanation  has  been  found  for 
future  efforts  should  be  concentrated  o 
discrepancy  between  theory  a 

The  steady  state  mean  vorticity  v 
following  three  ways: 


the  disagreement  a 
n resolving  the 


(a)  1 


e computations  based  o 


f rotation  of  th 
mean  vorticity 


(b)  by 


computations 


on  the  relationship 
t components,  and 


?n  vorticity  and  the  mean 


mill  and  the  computations  based  on  the  relationship  between 
vorticity  and  angular  velocity. 

The  results  from  these  three  approaches  were  in  general 


The  floating  disks  used  to  measure  mean  currents  did 
not  experience  any  rotation  outside  the  surf  zone;  however, 
rotation  started  as  soon  as  the  disks  entered  the  surf  zone. 
This  observation  agrees  with  the  theoretical  considerations 
that  mean  vorticity  is  only  generated  as  a result  of  energy 


the  mill  used  to  "measure"  vorticity 
are  not  completely  understood  in  spite  of  the  fact  that  the 
results  inferred  are  reasonable  when  compared  with  an 

Further  studies  are  needed  for  a better  understanding  of  the 
mechanics  of  the  mill  under  the  effect  of  the  wave  field. 

Future  studies  should  include  testing  with  floating 
mills  and  "fixed"  mills.  Also,  a larger  scale  model  would 
be  desirable,  as  a larger  surf  zone  would  allow  a more 
detailed  and  smoother  hydraulic  regime;  consequently,  better 
measurements  would  be  obtained.  The  mill  should  occupy  a 
larger  portion  of  the  total  water  depth  but  avoid  contact 
with  the  bottom  during  the  passage  of  the  trough.  Also, 
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measurements  of  the  mean  water  level  are  suggested  because 
currents  inside  the  surf  zone  are  driven  by  water  slopes. 
Finally,  an  automated  way  to  analyze  the  data  should  be 
developed. 

The  primary  purpose  of  this  work  was  not  to  study 
sediment  transport;  however,  the  following  comments  can  be 
made  based  on  hydrodynamic  consideration. 

(a>  As  a result  of  the  refraction  process,  convergence 


over  the  crest  and  later  over  the  trough  inducing  a set-up 
responsible  for  the  generation  of  currents  over  the  crest  of 
the  bar  to  shore  and  back  through  the  trough.  This  circula- 
tion associated  with  the  bar  might  be  interpreted  as  a self- 
maintenance mechanism  for  these  sand  features  as  indicated 
by  Niedoroda  { 1972)  after  a qualitative  analysis  of  natural 

with  Niedoroda's  observations,  for  those  bars  confined  to 

(bl  Since  these  circulations  have  a scale  proportional 
to  the  width  of  the  surf  zone,  no  conclusions  about  large 
scale  bars  (also  called  "sand  ridges"  which  have  a length  50 
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Table  A-l.  Wave  Height  Distribution  in  cm;  Wave 
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Table  A-3.  Wave  Height  Distribution  in  cm,  Wave 
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Table  B-5.  Angle  of  Rotation  8m  Versus  Time;  Run  Ho.  042; 
Depth  = ^3.2  cm?  Period  = 1.5  s;  Location 
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Table  B-8 . Angle  of  Rotation  9m  Versus  Time;  Run  No.  49; 

Depth  = 4.1  cm;  Period  = 1.0  s;  Location 
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Rotation  9m  Versus  Time;  Run  No. 
3.5  cm;  Period  = 1.0  s;  Location 
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Rotation  6m  Versus  Time;  Run  No.  56; 
3.7  cm;  Period  = i.O  s;  Location 
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Angle  of  Rotation  0m  Versus  Time;  Run  No.  57; 
Depth  = 4.7  cm;  Period  = 1.0  s;  Location 
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Table  B-26.  Angle  of  Rotation  8m  Versus  Time!  Run  No.  162: 
Depth  - 3.3  cm;  Period  - 1.5  s;  Location 
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sir  Is 


A.-,  ? , , f i 


K. 


-J- 


S*  A ST  A Sr  A Sr  A 


158 


TIME  (s) 

.?  i . f ',a 


>r 


S’ -60 


V 


-120 


159 


Rotation  9m  Versus  Time;  Run  No.  022; 
3.1  cm;  Period  = 1.5  s ; Location 


TIME  (s) 

(Aji-i  ? . i 6 , i ? , , '2,  | ■,$,  ,'8 

\ 

-30- 

V # 

S, -SC- 

OT 

-90- 


-I20l 


160 


;.r  a sr  ,.s,  sr  (deg)  sr  ..a, 


TIME  (s) 


Table  B-33.  Angle  of  Rotation  0m  Versus  Time;  Run  No.  160; 
Depth  *o4-1  cm;  Period  = 1-5  s;  Location 


TIME  (s) 

(ft. , i ? . i ? , , ? , . 'F  , , ',5  , ,'S 


-BO- 

'S 

f -SC- 

-90- 

L 

J 


162 

Table  B-34.  Angle  of  Rotation  6m  Versus  Time;  Run  No.  002- 
Depth  = 3.1  cm;  Period  = 1.5  s;  Location 


TIME  (s) 
o?  i i y , . f , , ? . . 

-30-  • 

5,-60  - 

5r 

-50- 


-I20l 


-6CL 


Angle  of  Rotation  8„  Versus  Time;  Run  No.  156; 
Depth  = ^3.7  cm*  Period  = 1.5  s;  Location 


dr^r-r-?- 


TIME  (s) 

9 » 


-30- 


-90 


-120 


Angle  of  Rotation  em  Versus  Time;  Run  No.  200; 
Depth  = 2.9  cm;  Period  = 1.9  s:  Location 


(deg) 


(deg) 


time  em 
Is)  (deg) 


(deg) 


(deg)  (s)  (deg) 


TIME  (s) 

ft-.?.  .9.  ,'P 


-30 


g.-60- 

-90- 


-120 


167 


sr  a,  sr  a sr  A Hr  A, 


3°r 


o^A-.  i i i j i ; i i i i i i i i 
3 .6  . 9 12  15  18 

TIME(s) 


-23- 


-80- 


-90 


Table  B-40.  Angle  of  Rotation  9m  Versus  Time;  Run  No.  208; 
Depth  = 3.3  cm;  Period  = 1.9  s;  Location 


-eoL 


169 


60r 


30- 


12  15  18 

TIME(s) 


-30- 


I Location 


y Floating 
k Trajectories;  Run  No.  41/01;  Wave 
iod  *ls;  Wavemaker  Starting  Time  = 


motion  becomes  perceptible. 


• ^^Initiol  F 


l. 


n Current  M 
k Trajectories?  Run  No.  41/02; 
iod  = 1 s;  Wavemaker  Starting  T 


Table  C-3.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories?  Run  No.  41/03?  Wave 
Period  = 1 s?  Wavemaker  Starting  Time  = 10.5  s 


U 


175 


n Current  Measurements  by  Floating 
k Trajectories!  Run  No.  42/01;  Nave 
' r Starting  Time  - 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


(m) 

* Trough 

L, 

1 

*- Breaker  Line 

" u ^Jlnitjal  Position 

, . ft, 

/Crest 

176 


n Current  Measurements  by  Floating 
k Trajectories  * Run  No.  42/02;  Wave 
iod  = 1 s;  Wavemaker  Starting  Time  = 


(m) 


‘Trough 


L, 


,-<t  | -Ores I | 

5 x(m) 


Position 


177 


n Current  Measurements  by  Floating 
k Trajectories*  Run  No.  42/03;  Wave 
iod  =ls;  Wavemaker  Starting  Time  - 


1 Trough  » 

—Breaker  Line 


*(m) 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories,-  Run  No.  43/01;  Wave 
Period  - 1 s;  Wavemaker  Starting  Time  = 16.2  s 


x(m) 


179 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  43/02;  Wave 
iod  = 1 s;  Wavemaker  Starting  Time  = 


Breaker  Line 
- — Initial  Fbsition 


-Crest 


n Current  Measurements  by  Float 
k Trajectories;  Run  No.  43/03; 
iod  = 1 s;  Wavemaker  Starting  T 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


' N~-  Initial  Position 


u 


ft  I /-Crest 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  43/04;  Have 
iod  = 1 s;  Wavemaker  Starting  Time  = 


^Trough  J 

— Breoker  Line 


— —Initial  Position 


ts  by  Floating 
No.  43/05;  Wave 
Starting  Time  = 16.2  s 


(a)  (m)  (m)  (cm/s)  (cm/s) 


motion  becomes  perceptible. 


Position 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  43/06;  Wave 
iod  = 1 s;  Wavemaker  Starting  Time  - 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


Table  C-13.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  44/01;  Wave 
Period  = 1 s;  Wavemaker  Starting  Time  = 29.5  s 


• ^Trough  j 
16  _ Iniliol  Ftsilicn  ' -i 


f Crest 


n Current  Measurements  by  Floating 
It  Trajectories;  Run  No.  44/02;  Wave 
iod  = 1 s;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


- — Initial  Position 


186 


s by  Floating 
s;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


'osition 

ft,  roa 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  44/05;  Wave 
Period  = 1 s;  Wavemaker  Starting  Time  ■ 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  44/06;  Have 
Period  = 1 s;  Wavemaker  Starting  Time  1 29.5  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


LJ 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  45/01;  Wave 
iod  -Is;  Wavemaker  Starting  Time  - 


£ Trough 

V 

I— , 

• • Initial  Position 

Breaker  Line 

i i ti. 

c Crest 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  45/02;  Wave 
Period  » 1 s;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


192 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  45/03;  Nave 
Period  = 1 s;  Wavemaker  Starting  Time  = 


£ Trough  v I 

•“-Initial  Position  [_ 


ft  i f Crest  I 


193 


ts  by  Floating 
Starting  Time  = 15.4  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


■ — Initial  Position 


-Crest 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  45/05;  Have 
Period  » 1 s;  Wavemaker  Starting  Time  = 15.4  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


^Trough 


^Initiol  Position 


ft  i fOest 


195 


ts  by  Floating 
Starting  Time  = 15.4  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


— Breoker  Line 


L. 


^ ^Jnit  iol  Position 
i*  \ i ft  i .'•Crest 


ts  by  Floating 
Starting  Time  = 15.4  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


u 


il  Position 

r<t  i r Crest  j 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  45/08;  Have 
Period  = 1 s;  Wavemaker  Starting  Time  ■ 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


^Trough 

L„ 

V u"* 

''X  Initial  Position 

- • • r"  V i d±_ 

rOest 

x(m) 


» 


Table  C-2B. 


n Current  Measurements  by  Floating 
k Trajectories;  Bun  No.  30/02;  Have 
iod  = 1.5  s;  Havemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


n Current  Measurements  by  Floating 
k Trajectories?  Run  No.  30/03;  Wave 
iod  = 1.5  s;  Wavemaker  Starting 


in  the  table  indicates  the  time  when  the 
perceptible. 


wris. 


202 


Mean  Current  Measurements  by  Floating  Disk 
Trajectories;  Hun  No.  30/05;  Wave  Period(s), 

' ‘ r Starting  Time(s),  23.2  seconds 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


1 * Trough 

u 

•-•^—Initial  Position 

'^Breaker  Line 

i 1 1 clt_ 

f Crest 

n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  30/06;  Have 
iod  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  30/07;  wave 
iod  »ls;  Wavemaker  Starting  Time  = 


.'WiTiS. 


207 


u 


ft,  fCresi 


S'  i« 

(m)  (cm/s)  (cm/s) 

Sill  !:i! 

K:K  |:]| 

!:!! 

S:!i  -3:!  :J:i 

!:Ji  :S 

1 1 |j 

|:|  |:j  l°l 

!!:i;  j|!i 

II  I 

II  1 

II  1 

!:il  :ll:i  :i:i 

!:S  -ll:i  !:S 

S:S  :l:l  !S 

1 i Jil 

1 1 :i 

!:!I  S:?  !:5 

Continued 


i — fn  fC:r:! 


.US'S, 


5:!  !:! 

-I  1 

-!:!  ":! 

i:l  J;: 

:l:l 

SJ  -1:1 

:5:S  :l:l 

1;! 

‘5:!  1:1 

S:i 

1:S  -3:! 

:l 

-3:S  !:S 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  31/06;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  31/07;  Wave 
iod  * 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


by  Floating 
. 32/Olt  Wave 
Starting 


Hi  ffr 


by  Floating 
Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


, £ Trough  j 

'Initiol  Position  I 


by  Floating 
. 32/03!  Wave 
Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


* Trough  \ 

- Inilicl  Rjsilion  T 


ft  i f Crest 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories?  Run  No.  32/04;  Have 
Period  = 1.5  s;  Wavemaker  Starting 


motion  becomes  perceptible. 


n Current  Measurements  by  Floatin 
k Trajectories;  Run  No.  32/05;  Wa 
iod  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible 


•-Trough 

L 

*•* — Iniliol  Position 

• • ' Breaker  Line 

/■Crest 

w* 


K:S 


ill 


-;S:I 


-1:1 


:iii 


Is 

!:! 


“:S 


laSl:  1 


Disk  Trajectories;  Run  No.  32/07;  Wa 
Period  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


'Ll 


^Breaker  Line 
f—  Initiol  Posilion 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  32/08;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


•* — Initial  Position 


225 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  33/01;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


v 

L_ 

'initial  Position 

'•^r-Breaker  Line 

, n ■ ft. 

/■Crest 

Disk  Trajectories; 


by  Floating 
Starting 


Time  x y 0 V 

(5)  (m)  (m)  (cm/s)  (cm/s) 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


/ 

iniliol  Position 


L 


227 


k Trajectories;  Run  No.  33/03;  Ha 
iod  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


\ • £Trough 

L 

\ *•  *^v“lniliol  Position 

1 5 i 

r Crest 

n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  33/04;  Wave 
iod  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


1 Trough 
ol  Position 


Breaker  Line 


°2 


; Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


^Trough 
— Initiol  Position 


230 


by  Floating 
. 33/06;  Wave 
Starting 


(s.)  (m)  (m)  (cm/s)  (cm/s) 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


\ereoKer  Une 

u 

. : —'"i 

ol  Position 

ft, 

/Crest 

n Current  Measurements  by  Floating 
It  Trajectories;  Run  No.  33/07;  Wave 
iod  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


Initiol  Posilion 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories?  Run  Mo.  33/08;  Wave 
Period  = 1.5  s;  wavemaker  Starting 


motion  becomes  perceptible. 


£Trouqh  Y 

— Breaker 


^-Initial  Position 

, ffci  /-Crest 


xlm) 


233 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  33/09;  Have 
Period  = 1.5  s;  Wavemaker  Starting 


The  £irst  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


u 


|*LFPS"HHs'lE?- 


Continued 


‘■Trough 


Initial  Position' 
v "^Breaker  Line 


u 


EJfspS-SS* 


n Current  Measurements  by  Floating 
k Trajectories,-  Run  No.  34/03i  Wave 
iod  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


n Current  Measurements  by  Floating 
k Trajectories!  Run  No.  34/04;  Wave 
iod  = 1.5  s;  Wavemaker  Starting 


Table  C-62.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  34/05;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


^Trough 


L 


Initial  Position 


ft.  f Crest 


Table  C-63.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  34/06;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


“ gigSS-2^ 


by  Floating 
. 34/08;  Wave 
Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


(A) 


o; 


'f'Breaker  Line 

'T"',h  [_u 

\ . . 

• *i ^ fti  r Crest 

n Current  Measurements  by  Floatin 
k Trajectories!  Run  No.  34/09;  Wa 
iod  = 1.5  si  Wavemaker  Starting 


^Breaker  Line 

£ Trough 

u 

Initial  Position 

ft, 

f Crest 

3 4 • 

*(m) 

n Current  Measurements  by  Floating 
It  Trajectories;  Run  No.  35/01;  Wave 
iod  = 1.5  s;  Wavemaker  starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


/' 


Initial  Position 
^Breaker  Line 


t_ 


Table  C-60 . Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  35/02;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  35/03;  Wave 
iod  = 1.5  s;  Wave maker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


H1Sv“ 


rr'W.sra. 


(« 

(cm/s) 

(cm/s) 

j£jj 

!4j 

‘if! 

4:1 

:Sil 

SS4S 

!:S 

44 

**ji 

Sis 

;iit 

lil 

if 

sis 

1:1; 

si 

is 

}";*« 

!:!! 

4?:S 

ul:°A 

S:S 

-45:1 

The. first 

time  in  the  t 

es  the  time  wh 

en  the 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  35/OS;  Wave 
Period  = 1.5  s;  Wavemaker  Starting 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


\ 

IT”'“ 

L 

Initial  Position 

3 

• Breaker  Line 

i i 

f Crest 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  35/06;  Have 
Period  = 1.5  s;  Wavemaker  Starting 


Position 


:-73. 


gffirHS-HSr- 


S:!! 


IS 


lia 

;| 


i:H 


1 


3:2 


sr 


1:51 

I: 


sis 


!:5! 


*1 


!;y 

1 

sis 


t: 

: 

I; 

3:i 

Si 

sil 

a 


252 


Continued 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


.bw!ra. 


is 

HI 
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A 

it 
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HI 
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HI 
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HI 
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Hi 

S3 
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Hi 

11 
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table  indicat 

es  the  time  • 
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Continued 


“ Inf  5=ss-sa?« 


(m)  Icm/Bl  (cm/s| 


Sr 


K-ur1; 


SJSS-..,.. 


!£  SiB 

! I 

is  Sis 

II  if 

1:1!  i:ll 


Si 


:l:S 

Si 


V, 

Si 


Si 


(m) 

Breaker  line 

^Trough  | 

i Ni  i 

__i ft.  fc”' 

Table  C-79.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  01/07;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 20.9  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


£ Trough  v 


,-ti  .'■Cresi  1 


sssssrs's.wjr;:;. 


r ' \ i i ft  i [Crest 

3 4 5 6 x(m) 


ktitikdzz 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  02/03;  Wave 
Period  = 1.9  s;  Wavemaker  Starting  Time 


<s)  (m)  (m)  (cm/s)  (cm/s) 


The  first  time  in  the  table  indicates  the  time  when  the 
• motion  becomes  perceptible. 


Breoker  Line 

! Trough  y 

L-n 

( . 

• • • •—Initial  fiisition 

i ^ i 

i rt  ■ f Crest 

x(m) 


n Current  Measurements  by  Floating 
k Trajectories?  Run  No.  02/04;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


j 

^Trough 

u 

'* — Breoker  Line 

• • ——Initial  Position 

• i ’ ' , 

l ft, 

[ Crest 

0I ll  1 I I .-Cres) 

2 3 4 5 6 x(m) 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  02/05;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  02/06;  Wave 

r Starting  Time  = 7 


V* — Breaker  Line 

sTrough 

L 

.Initial  Position 
-■  i 

/Crest 

n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  03/01;  Wave 
iod  - 1.9;  Wavemaker  Starting  Time  = 


motion  becomes  perceptible.  time  Wh6n  th6 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  03/02;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  - 


u 


267 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  03/03;  Have 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


L„ 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  03/04;  Have 
Period  = 1.9;  wavemaker  Starting  Time  = 11.0  s 


motion  becomes  perceptible. 


V* — Breaker  Line 


Position 


n Current  Measurements  by  Floating 
k Trajectories?  Run  No.  03/05;  Wave 
' ' ’ * " " ‘ r Starting  Time  = 


motion  becomes  perceptible. 


8reoker  Line 


^Trough 


it  i [ Pest 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  03/06;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  - 


1 t* — Breaker  Line 


■Initial  Position 


271 


y Floating 
k Trajectories;  Run  No.  04/01;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


u 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  04/02;  Have 
r Starting  Time  * 


motion  becomes  perceptible. 


i— Breaker  Line 
1 

■ .*•— Initial  Position 


irougn  | 

ffci  f Crest 


by  Floating 
oriest  Run  No.  04/03;  Wave 
; Wavemaker  Starting  Time  * 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


— Initiol  Position 


u 


ft  I f Crest 


Table  C-95.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  04/04;  wave 
Period  = 1.9;  Wavemaker  Starting  Time  - 12.9  s 
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n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  04/05;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


r*—  Breoker  Line 


u 


°r 
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Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  04/06;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


— Breaker  Line 


L. 


'•v.V 


■Initial  Position 


ft.  i --Crest 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  04/07;  Have 
Period  = 1.9;  Wavemaker  Starting  Time  = 12.9  s 


(••—Breaker  Line 


u 


ft  I fCresl 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  05/01;  Wave 
iod  - 1.9;  Wavemaker  Starting  Time  = 


<s)  (m)  (m)  (cm/s)  (cm/sl 


motion  becomes  perceptible. 


(m) 

'r — Breaker  Line 

£Tiouqh  | 

^ ^Initial  Position 

i ft;  fCrest 
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Table  C-100. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  Mo.  05/02;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  - 10.5  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


— Breaker  Line 

•* — Initial  Position 


^Trough 


u 


6 x(m) 


k Trajectories;  Run  No.  05/03;  Wave 
iod  - 1.9;  Wavemaker  Starting  Time  = 


. J*—  Breaker  Line 

£Trough  j ^ 

'O1- Initial  Position 

i ft  i f Crest 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  05/04;  Have 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


h — Breoker  Line 

**t  Initial  Position 


■ fCresl 


:-103.  Continued 


k Trajectories;  Run  No.  05/06;  Have 
iod  = 1.9;  Wavemaker  Starting  Time  = 


u 


ft  I fCresI 
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Mean  Current  Measurements  by  Floating 
Disk  Trajectories!  Run  No.  06/01;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


1 

1 Trough 

L 

• t—  Breoker  Line 

" • ^--.^^Initiol  Position 

ft, 

r Crest 

n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  06/02;  Wave 
iod  » 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


4T 


Breaker  Line 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  06/03;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


ts  by  Floating 
Ho.  06/04;  Wave 
Starting  Time  = 19.1  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


Table  C-109.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  06/05;  Mave 
Period  = 1.9;  Wavemaker  Starting  Time  = 19.1  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


^Trough  \ 

- Breaker  Line  L .u 


"•^Initial  rositic 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  06/06;  Have 
iod  = 1.9;  Wavemaker  Starting  Time  = 


C-lll.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  06/07;  Have 
Period  = 1.9;  Wavemaker  Starting  Time  = 19.1  s 


(si  (">)  (m)  (cm/s)  (cm/s) 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


Position-.  A 


u 

ft  i f crest 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  06/08;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  3 


w — Breaker  Line 

* Trough  v. 

1 -u 

Position^  i ''i  i 

! T Crest 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  06/09;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


[Posmony  rCres. 

2 . . ••  • "3  4 5 6 x(m) 


C-114.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  Mo.  07/01;  Have 
Period  = 1.9;  Wavemaker  Starting  Time  = 25.8  s 


motion  becomes  perceptible. 


V*—  Breoker  Line 


L 


rt  i /Crest 


Table  C-115.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  07/02 ; Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 25.8  s 


IT 

w A 

(cm/s) 

(cm/s) 

jo  :jj 

3!l7  o!oO 

'-IX 

2^68  -0.24 

2.53  -0.27 

-5 ! 3 

time  in  the  table  indicates 
comes  perceptible. 

* 

Breoker  Line 

£Trough 

-L, 

. .*  '^/‘-■Initial  Position 

/ i . i 

ft. 

c Crest 

3 4 5 

6 

x(m) 

n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  07/03;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


(m) 

'''"-—Breaker  Line  tTrouqh 

L, 

Initial  Position 

/ . - . ft. 

fCrest 

1 1 ' t | i /-Crest 

2 3 4 5 6 x(m) 
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n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  07/04;  Wave 
iod  * 1.9;  Wavemaker  Starting  Time  = 


'"•Mr — Breaker  Line  tTrou9h 

L, 

•*  * * '^'^’‘-'Initiol  Position 

! . , , ft. 

c Crest 

~'z 


n Current  Measurements  by  Float 
k Trajectories;  Run  No.  07/05; 
iod  = 1.9;  Wavemaker  Starting  T 


"'''^—Breaker  Line  ^Trough 


Table  C-119.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  07/06;  Wave 
Period  = 1.9;  wavemaker  Starting  Time  = 25. B s 


' ——Breaker  Line 


L, 


C ti  fCresl 


It  Trajectories;  Run  No.  07/07;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


301 


RunnNl.b08/S“tWave 


1 

1 


303 


Table  c-122 . Continued 
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»r 


iTsr; 


!S:!! 


i:S 

:1 

!:?! 

I 


H 


"/*"  ;Troliqh 


u 


ft  I fCresl 


Table  C-127. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  09/02;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 15.7  s 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  09/03;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


The  tirst  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


y j 7 " * ^Trough.  *" 

■m'  | Initial  Position  • 

V*  Breaker  Line 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories!  Run  No.  09/04;  Wave 
Period  = 1.9'j  Wavemaker  Starting  Time  * 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


Table  C-130. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  09/05;  Wave 
Period  a 1.9;  Wavemaker  Starting  Time  = 15.7  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


n Cuccent  Measurements  by  Floating 
k Trajectories;  Run  No.  09/06;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 1 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


1 Initial  Position—" 

^Trough  * 

, 

: L- 

p — Breaker  Line 

v 

, r t t /Crest 

n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  09/07;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = IS. 7 


n Current  Measurements  by  Floating 
It  Trajectories;  Run  No.  10/01;  Wave 
r starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


^ Tough  \ 


ft  I fC.es!  | 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  10/02;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  ■ 


motion  becomes  perceptible. 


•—Initial  fijsition  £ Trough 


r— Breaker  Line 


L. 


ft  i f Cfesl 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  10/03;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


•* — Initial  Rjsition  ^Trough 

L. 

■ 1 

• • . , i ft. 

/Crest 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  10/04;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


(£> 

Initial  Position  ^Trough 

L 

V“ — Breaker  Line 

\ 

■ » £ 1 i_ 

3X9 


Table  C-137. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  10/05;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 8.1  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


— Initial  Position  ^ Trough 


L 


Table  C-138. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  11/01;  Wave 
Period  = 1.9:  Wavemaker  Starting  Time  = 15.5  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  11/02;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  = 


•*—  Initial  Position  1 Trough 


L» 


xlm) 
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n Current  Measurements  by  Float 
k Trajectories;  Run  No.  11/03; 
iod  = 1.9;  Wavemaker  Starting  T 


“Trough  » 

L, 

ft  i fCresi 


Breaker  Line 


C-141.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  11/04;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 15.5  s 


u 


Table  C-142. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  12/01;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 11.0  s 


x! 


n Current  Measurements  by  Floating 
k Trajectories;  Run  No.  12/02;  Wave 
iod  = 1.9;  Wavemaker  Starting  Time  ■ 1 


• Breaker  Line 


Table  C-144.  Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  12/03;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 11.0  s 


I* — Breoker  Line 

^Trough 

L. 

y •* — Initial  Position 

ft, 
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Table  C-145. 


Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  12/04;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 11.0  s 


•s*  <m>  (s'!  (cm/s)  (cm/s) 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


t* — Breoker  Line 

t Trough 

Lt 

^ •*—  Initial  Position 

. ft, 

/■Crest 

3 4 

5 6 

<(m  ] 

Mean  Current  Measurements  by  Floating 
Disk  Trajectories;  Run  No.  12/05;  Wave 
Period  = 1.9;  Wavemaker  Starting  Time  = 11.0  s 


The  first  time  in  the  table  indicates  the  time  when  the 
motion  becomes  perceptible. 


fCres! 


Dye  Releases 


147.  Mean  Current  Measurements  by 


Mean  Velocity 


(cm/s) 
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Mean  Current  Measurements  by  Dye  Releases; 


(cm/s) 


(cm/s) 


APPENDIX  D 

"VORTICITY  METER":  THEORETICAL 

EXPERIMENTAL  CONSIDERATIONS 


iglSpfe, 


“Mill”-, 
Meter  Rod-^ 

f 

— Support 

LATERAL  VIEW 
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Position 


Mill  in  Its  Natural  (No  Extra  Load) 
Position  of  Equilibrium;  P Represents 
the  Net  Height  of  the  Mill 


mill 


5! 

S 

S! 


!:S! 

i:IS 


&(degre 


Angle  of  Rotation,  $£,  Versus  Time  as 
the  Kill  Recovers  Its  Position  of 
Equilibrium 


341 


Estimation  of  the  Angle  of  Rotation  of  the 


A limited  number  of  tests  was  carried  out  in  order  to 
calculate  the  rotation  of  the  mill  under  different  wave 
conditions.  The  experiments  were  performed  at  two  different 
depths  of  measurement,  wave  heights  and  periods.  The  angle 
of  rotation  was  recorded  on  video  tape.  As  an  example.  Test 
No.  2 was  performed  (location  of  the  mill)  15  cm  from  the 
water  surface,  the  steady  state  wave  height  was  0.6  cm,  the 
wave  period  1.3  s and  no  continuous  rotation  was  observed. 

The  angle  of  rotation  versus  time  for  each  of  these 
experiments  are  presented  in  Tables  D-6,  D-7  and  D-8. 

Figure  D-6  shows  the  angle  of  rotation  of  the  mill  versus 


the  first  180  degrees  of  rotation  the  natural  torque  of  the 
mill  opposes  the  wave-induced  torque  and  for  the  next  160 
degrees  of  rotation  it  works  in  favor  of  the  wave-induce 
torque  increasing  the  angular  velocity,  represented  by  the 
slope  of  the  curves. 


Table  D-5.  Angle  of  Rotation  of  the  Mill,  0m, 


Table  D-6.  Angle  of  Rotation  of  the  Mill,  8n 


of  Rotation  (0) 
(degrees) 


Table  D-7.  Angle  of  Rotation  of  the  Mill,  0m, 
Versus  Time,  Test  No.  10 


of  Rotation 
(degrees) 


(0) 


Angle  of  Rotation  of  the  Mill,  0m, 
Versus  Time,  Test  No.  13 


Time  Angle  of  Rotation 

Is)  (degrees) 


(0) 


rotation  was  not  observed,  the  limited  rotation  was  in 
static  equilibrium  with  the  wave-induced  kinematics.  In 
order  to  compute  the  wave-induced  torque  for  those  cases, 
the  known  natural  moment  of  the  mill  was  used.  As  an 
example,  consider  Test  No.  1.  The  mill  rotated  65  degrees 
before  reaching  equilibrium  with  the  wave  field,  therefore 
the  averaged  wave-induced  torque  is  • 


350- 


300- 


250- 


100  - 
50- 


• Test  No.  3 
° Test  No.  13 


0 


10  20  30 

TIME  (s) 


Angie  of  Rotation  of  the  Mill,  0,  Versus 
Time  under  the  Effect  of  a Wave  Field; 
Test  No.  3;  Wave  Height  = 5.5  cm; 

Period  = 1.0  s;  Depth  = 48.0  cm;  Depth 
of  Measure  = -15.0  cm;  Test  No.  13;  Wave 
Height  = 12.0  cm;  Period  = 1.1  s; 

Depth  = 48.0  cm;  Depth  of  Measure  =*  -36.0 


"1 


the  mill 
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Syi 


(Fu+,  Fu _)  represent  the  x-components  o£  the  drag  force  on 
the  upper  and  lower  vanes , respectively,  and  (Fw+,  Fw_) 
represent  the  z-components  of  the  drag  force  on  the  left  and 
right  side  vanes,  respectively.  The  differential  form  of 


where  CD  is  the  drag  coefficient,  p the  water  density,  b the 
width  of  the  vane  and  dr  represents  a differntial  length 
along  the  vane,  (u+,  u_)  the  local  horizontal  velocities  on 

cal  velocities  on  the  left  and  the  right  vanes  under  the 
assumption  that  the  mill  does  not  rotate.  The  instantaneous 
torque  (Tq) , in  its  differential  form,  is  defined 


velocities 


expressed 


where  (u  , wc)  represent  the  velocity  component  at  s = s 
Substituting  the  above  in  Equation  (D-9) , integrating  in 
and  time  averaging  over  one  wave  period  the  net  torque 


- <ue 


+ (we 


(D 


Before  proceeding  further,  consider  a 
example,  the  instantaneous  applied  torque 
the  wave.  The  integrand  of  Equation  ID-11 


simplified 


instantaneous  torque 


Based  on  linear  wave  theory  the  wave  orbital  velocities  are 


and  the  gradients 


Substituting  in  Equation  (D— 12)  and  considering  Equation 


where  the  dispersion  relationship 


5-18) 


has  been  considered.  Figure  D-9  displays  the  nondimen- 
sional  wave-applied  torque  under  the  crest  of  the  wave 
(assuming  that  the  mill  does  not  rotate)  for  different 
values  of  h/L. 

Consider  now  the  net  applied  torque  (averaged  over 
wave  period) . From  Equations  (D— 1 4 ) and  ( D— 1 6 ) the 
integrand  in  Equation  (D-ll)  results  proportional  to 


|cc 


- ot)  | 


which  is  identically  equal  to  zero  when  averaged  over  one 
wave  period.  Therefore,  linear  waves  do  not  contribute  to 
the  net  torque  and  only  induce  an  oscillatory  torque  which 


Consider  now,  the  effect  of  the  nonlinearities  on  the 
net  torque  assuming  a wave  field  represented  by  a second 
order  Stokes  wave.  Under  this  wave  the  orbital  velocities 


uc  “ °1  ch  lk  so>  cos  ,kx  ‘ 041 

♦ U 2 ch  (2k  sD)  cos  2(kx  - ot) 


°2  = W2 


TT  ■ °lkshlk  so>  “s  «*  - #t» 

+ °2  2k  sh  (2k  so>  cos  2(k*  - ot> 
T3T  = H1  k sh  (k  se)  cos  (kx  - .t) 

+ W2  2k  sh  (k  sQ)  cos  (kx  - at) 


As  a numerical  example  consider 
D-4.  The  wave  characteristics  v 


and  the  other  parameters  used  in  the  calculations 

drag  coefficient  (CD)  : 1.0 

water  density  (c)  : 1.0  gram/cm 

length  of  the  vane  (R)  : 4.3  cm 

The  resulting  maximum  velocities  were 


the  strength  of  the  nonlinearities  can  be  expressed 


velocity  component  amplitude 


Figure  D-10  shows  the  velocity  components  as  a function  of 
time.  To  calculate  the  Integral,  the  fourth  order  Simpson's 
rule  was  applied.  The  results  are  summarized  below 


net  torque  on  the  lower  vane 
net  torque  on  the  left  vane 


The  net  torque  obtained  is  two  orders  of  magnitude 
smaller  than  the  minimum  torque  needed  to  induce  a 
continuous  rotation.  As  a consequence,  based  on  this 
analysis  it  is  not  possible  to  explain  the  rotation  of  the 

gradients  in  the  wave  height,  the  following  considerations 
are  necessary:  before  a broken  wave  arrives  to  the  place 
where  the  mill  is  located,  the  wave  height  gradient  (in  the 
y-direction)  is  negative,  so  larger  wave  heights  are  on  the 
crest  of  the  "bar"  due  to  refraction,  inducing  a clockwise 
rotation  (negative  vorticity) . When  waves  break  the  wave 
height  gradient  becomes  positive  because  the  wave  height, 
inside  the  breaker  zone,  is  proportional  to  the  local  water 
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depth,  consequently  induces  a positive  rotation. 

Measurements  of  the  angle  of  rotation  of  the  mill  indicate 
negative  values.  Therefore,  it  is  clear  that  the  wave- 
induced  vorticity  overcomes  the  extraneous  effect  introduced 
by  the  gradient  of  the  wave  height. 
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